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Abstract 

A general class of nonlocal cosmological models is considered. A new method 
for solving nonlocal Friedmann equations is proposed, and solutions of the Fried- 
mann equations with nonlocal operator are presented. The cosmological properties 
of these solutions are discussed. Especially indicated is p-adic cosmological model 
in which we have obtained nonsingular bouncing solution and string field theory 
£T) . tachyon model in which we have obtained full solution of nonlocal Friedmann equa- 

tions with w = — 1 at large times. The possibility of obtaining realistic value of 
cosmological constant from nonlocal cosmological models is also discussed. 

in 

i> ■ 1 Introduction 

o 

■—^ In recent works there appears interest in nonlocal cosmological models derived from 
String Field Theory in connection with problem of describing cosmological inflation or 
accelerating expansion of the Universe. 

Modern cosmological data indicates that expansion of the Universe is accelerating. 
It may be owing to a component of the Universe with negative pressure, Dark Energy. 
Recent results of WMAP [1] together with la supernovae data give us the following range 
for the dark energy state parameter w = — 0.97^QQg. In the assumption that the value 
of the state parameter changes in the range > w > — 1 there exist different theoretical 
models for the dark energy (see reviews [2, 3] and references therein) - quintessence 
models [4], K-essence [5], quintom models [6], models constructed on the DBI action 
(see reviews [7] and references therein, [8]) and dilatonic models [9]. The case w = —1 
corresponds to models with cosmological constant, which are also possible. We also can 
find in the literature models related w < — 1 with modified gravity [10]. 

A nontrivial possibility is the case w < — 1, that means the violation of the null 
energy condition (NEC). One of the possibilities is to consider the phantom Universe 
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[11, 12]. There were proposed several coupled scalar-gravity models in which the null 
energy condition is violated, but as an effect most such models are unstable [13]. Recently 
it was proposed another possibility to consider dark energy description in the context 
of D-brane decay in cubic superstring field theory [11]. It was shown that D-brane 
decay can be interpreted at least at late times as a phantom model. Note that unlike 
phenomenological phantom models here phantom is an effective theory. Since string field 
theory by itself is a consistent quantizable theory this approach does not suffer from 
instabilities which are inevitable for phenomenological phantom models [11]. Such a 
nonlocal cosmological models were also considered in [14-18]. 

At the same time there have been a number of attempts to realize description of the 
early Universe via nonlocal cosmological models [18-21]. One example is p-adic inflation 
model [19] which is represented by nonlocal p-adic string theory coupled to gravity. For 
this model, a rolling inflationary solution was constructed and the interesting features 
were discussed and compared with cosmic microwave background (CMB) observations. 
The possibility of obtaining large nongaussian signatures in the CMB has also been 
considered in a general class of single field nonlocal hill-top inflation models [20] . Another 
example is investigation of the inflation near a maximum of the nonlocal potential when 
non-local derivative operators are included in the inflaton Lagrangian. It was found that 
higher-order derivative operators in the inflaton Lagrangian can support a prolonged 
phase of slow-roll inflation near a maximum of the potential [21]. 

It may happen that both early stage and contemporarily Universe can be explained 
with a single nonlocal model derived from string field theory. To be able to tell more 
about the early stage we need to consider perturbations in nonlocal models. It will be 
very interesting to compare obtained results in particular with those in [22-24]. To be 
able to do it the new technique for solving nonlinear Friedmann equation is required 
which will hopefully be provided in the present work where for the first time construction 
of solution of the full nonlocal Friedmann equations is presented. That is why in the 
present work we discuss numerical solution of Friedmann equations and their features in 
more detail. 

One of the popular examples of nonlocal models is the D-brane decay in the grav- 
itational background in the framework of string field theory which is described by the 
following tachyon action 

d'xV^ 1-fR + \4n g <j> + -0 2 — -$ 4 — t — A'J , (l) 

where g is the metric, m p is reduced Planck mass, <fi is tachyon field, $ = e kn <ft, k — |, 
□ g is usual Beltrami-Laplace operator, T is tension of the D3 brane, A' is the effective 
cosmological constant. 

Note that covariant string field theory as of today is constructed only in the flat back- 
ground. It was built only for some special types of non-flat backgrounds, in particular, in 
the anti-de sitter background [25] and in the plane wave background [26]. So the action 
([1]) at the special choice of V(<&) is a direct generalization for the approximated tachyon 
action to the case of arbitrary metric. 
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The action (pQ) in the Friedmann-Robertson- Walker (FRW) background 

ds 2 = —dt 2 + a (t)(dx\ + dx\ + dx\) 
leads to the following Friedmann equations 

3H 2 = — 2 S, (2) 

3H 2 + 2H = - At V, (3) 

which have the same form as the usual equations except for extra term in the expres- 
sion for the energy £ and pressure V which appear from the nonlocal interaction and 
which contain the Hubble function. This term makes these Friedmann equations more 
complicated even for numerical consideration. More precisely, 

£ = £ k + £ p + T + A' + £ nll + £ nl2 , (4) 

V v ' 

usual local energy expression 

V = £k - £ P - T - A' - £ n n + £ ni2 , (5) 

V v ' 

usual local pressure expression 

where 

£mi = k [ dp (e kpn ^ 3 ) (-n g e- kpa s$) , £ nl2 = -k f dp (de kpD ^ 3 ) {de~ kpD ^) . 
Jo Jo 

The field equation reads 

(£ 2 n 9 + l)e- n *$ = V'{$), 

here and below we will denote by / usual function derivative. The Beltrami-Laplace 
operator takes the form D g = —d 2 — 3H(t)d t + ^d 2 . 

In the case of cubic potential the action (jTJ) corresponds to the the tachyon in bosonic 
String Field Theory [27] for the lowest level in the level truncated scheme [28,29] with 
the induced twisted space-time. It is supposed that that we are dealing with D3-brane 
in the D26 space-time and the volume of the D22 compactified subspace is omitted. 

The case of quartic potential in the action (pQ) corresponds to inclusion of the metric 
in the tachyon action of the fermionic cubic string [30] in the approximation of a slowly 
varying auxiliary field [31]. Such an action appears in the level truncated scheme where 
only lowest scalar field (tachyon field) in the GSO(-) sector and correspondingly the 
lowest scalar field in the GSO(+) sector (this field does not have kinetic term and is 
considered as auxiliary field) are considered. Integration over the auxiliary field results 
in the quartic potential for the tachyon field. Note that appearance of the non-extremal 
Z)-brane in the framework of fermionic string field theory requires GSO(-) sector in the 
spectrum [32,33]. 

The purpose of this work is to consider a general class of nonlocal cosmological models 
of the form ([!]) and discuss which physical properties and consequences might be obtained 
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from such models. Among them we will be interested firstly in the classical solutions of the 
corresponding Friedmann equations which can be considered as a first approximation to 
the quantum solutions and might be useful for the study of ways to avoid the cosmological 
singularity problem [34]. 

We start with the description of some general class of nonlocal models. In section 
3 we will describe new numerical algorithm for solving nonlocal Friedmann equations. 
In section 4 we will present the numerical solutions of Friedmann equations for different 
parameters and will provide analytical explanation for some of their interesting features. 
In section 5 we will consider which cosmological solutions can be provided by nonlocal 
models. 



2 General Class of Nonlocal Cosmological Models 

In this paper we will consider the following nonlocal scalar field on the Z)-brane coupled 
to the gravity 



/ 



d 4 x\/—g 



M 2 M A ft 2 i \ 



(6) 



where g is the metric, D g = ■^=d ( i\f—gg ia 'dv- l M v is a Planck mass, M s is a character- 
istic string scale related with the string tension a' as M s = is a scalar field, 
#4 is a dimensionless four dimensional effective coupling constant related with the ten 
dimensional string coupling constant go and the compactification scale, A = ^-A' is an 
effective four dimensional cosmological constant, T is the brane tension, A; is a nonlocal 
coupling constant [31,33]. 

For further investigation let us rewrite our action in dimensionless space-time variables 



(7) 



S = J d^x^Tg {^R + ^0D 3 + ]f - Vi?) - A' - , 

where is a dimensionless scalar field, $ = e kUa (j) and m 2 — g^i- 
As a particular metric we take FRW one 

ds 2 = -dt 2 + a 2 {t) (dx\ + dx\ + dxf) . (8) 

We will consider spatially homogeneous configurations for which Beltrami-Laplace 
operator takes the form O g = —d 2 —3H(t)dt, for the convenience of numerical calculations 
let us introduce the following notation = d 2 + 3H(t)d t , which we will understand 
as generalization of second order time derivative, which contains extra term 3H(t)d t 
associated with gravitational background. 

Equation of motion for the space homogeneous configurations for the scalar field $ 
takes the form 

(-£ 2 2)h + l)e 2fcI) H$ = v' (<!>). (9) 
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The Friedmann equations have the following form 



3# 2 = ±e, 



m 2 



' 1 (10) 



3H 2 + 2H= V. 

For the case of nonlocal potentials the energy and the pressure have additional non- 
local terms S n n and S n i 2 

S = S k + S p + T + M + £ nll + S nl2 , (11) 

N v ' 

usual local energy expression 

V = £ k -£ p -T-M - S nll + S nl2 , (12) 



usual local pressure expression 



where 



£fc = |W, S p = -\<P 2 + V^), (13a) 
Sun = k f Q d P (e- kpV »V'(<f>j) (T> 2 n e k ^) , (13b) 

S nl2 = -kj dp (r^l/^)) (de kf)V ^ . (13c) 

As we can see the structure of Friedmann equations and of equation for scalar field 
in this model is rather complicated and their study is already a very interesting mathe- 
matical problem and can be considered as a separate mathematical investigation. 

To avoid calculation of e~ fcp:E, H term which is much harder to compute than e fcpI) H 
(k > 0) as computation of the former results in an ill-posed problem [35] we will use 
the following representation for nonlocal energy terms on the equation of motion for the 
scalar field 

£ n ii = k [ d P (VeX + l)e( 2 -^H$) [ple k ^) , (14a) 
S nl2 = -k jf * dp (d{-eT>l + l)e^ kV »^ (de k ^<b) . (14b) 

3 Method for Solution of Nonlocal Friedmann 
Equations 

3.1 Iterative Procedure for Solution Construction 

For numerical calculations let us rewrite our system in the following form 



1 



3H 2 = — S, (15a) 
m? p 



(-£ 2 D| + l)e 2fcI) H$ = V'($). (15c) 

In this paper we would like to point out which physical results we can extract from 
solutions of this system. First, we will find Hubble function and the scalar field from 
equations ( 115bl) and ( 115cl) . It happens that equation ( 115al) plays a role of energy conser- 
vation from which the effective cosmological constant is extracted and is unique for every 
field configuration. 

Integrating equation (115bl) we get the following system 

(-£ 2 D| + l)e 2fcI) H$ = v'(<$>), (16a) 

rc2 



H = / dr 



mj jo 



(16b) 



L^f -kj^dp (d(-£ 2 2)| + l) e ( 2 -^ fcI, H$) (de^W^ 

To find solutions of the system (fT6l) we construct the following iterative process 

$ = lim $ n , H = lim H m (17a) 

where iterations $ n and H n are in turn obtained as limits of sub-iterations 

= lim $ nm , if n+ i = lim H n>m , (17b) 

m— »oo m— >oo 

which are recursively defined as (m ^ 0) 

<Wi = + (17c) 

1 A , (9e *^ n , TO$n+1 )2_ (17d) 



H n , m +i = 5- / 

™p Jo 



kfdp (d{-evl n m + l)e (2 -" )fcI, H„, m<l , n+1 ) (9e^H-$„ +1 ) 

where initial iteration in m is taken as 

$„ i0 = $ n , F n>0 = F n . (17e) 

By T/ /_1 we denote a function inverse to V. Note how sub-iteration $„ )m depends only 
on H n , i.e. the whole step of sub-iterations for H. Note also that H n ^ m depends on $ n+ i, 
this is essentially a way to accelerate convergence. 

Intuitively this iteration procedure does the following. First fixing H n it finds $ n+ i us- 
ing iteration process (I17cl) which is a natural generalization of a process which is described 
in [36-38] and for which convergence is proved analytically. Then H n+ i is obtained by 
iterating f!17dl) where $ n+ i is fixed as computed in the previous step. Thus found H n+ i 
is used to find $ n+ 2 and so on. This kind of multi-step iterative procedure are known to 
provide rapid convergence. In our experiments this iterative procedure converges much 
more reliably as compared to more naive iterative procedure where one does not do 
intermediate iterations. 

In a similar fashion one can construct iterative procedures for systems with several 
scalar fields such as those studied in [35,39]. 
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3.2 Exponentiating by Solving Diffusion Equation 

Equations of motion as well as the iteration method discussed above require computation 
of action with operator e fcpI)H on a given function. Computationally this is a nontriv- 
ial problem by itself. We use the following method which beside giving mathematical 
definition for such an operator more importantly provides corresponding computational 
method. 

Result of the acting with operator e pl) H on a given function (p(t) is natural to define 
as a solution of the following diffusion partial differential equation with corresponding 
initial and boundary conditions 

= 2>h/(p,*), (18) 

f(0,t) = <p{t), (19) 
/(p,±oo) = <p{±oo), (20) 

where — oo < t < +oo, p ^ 0. Note that here evolution goes along p axis, while t plays a 
role of space. Solution of ([TBI is a rigorous mathematical definition of the operator e pI, H 
acting on function if 

ef"*<p(t) = f(p,t). 

Such solution exists and is unique for wide classes of H and if [40]. Since in standard 
settings such a diffusion problem is well defined in particular for H,ip e which 
seem to suffice applications which are considered here we will not discuss this topic any 
further. 



3.3 Fixing Potential and Other Parameters 

So far we did not specify exact shape of the potential for our system. In this paper we are 
interested in coupled to the gravity a tachyon field which appears as a lowest excitation in 
level truncation scheme for fermionic string field theory [31]. Thus all further discussion 
will be about quartic potential 

vm = V. (21) 

We would like to note though that from numerical perspective this potential can lead to 
difficulties as it results in solutions which are not different iable at zero (they behave as 

t l/3 

at in the vicinity of zero [41]). For the later reason we found it useful to study similar 
problem in potential 

\/ Q ($) = (l-a)i$ 2 + a i$ 4 , 0<a<l. 

As a — > 1 this potential tends to potential (1211) . The interesting property of this potential 
is that solution has finite derivative at zero and moreover there is a theorem which states 
convergence of our type of iterative procedure in case of Minkowski metric [42]. 
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Let us return to the potential f[2XT) - Equation (|15c[) is invariant under shifts 

$(t) + H(t)^H(t + t ), (22) 

and mirroring 

H (t)^-H(-t). (23) 

Thus it is natural to look for odd solutions. Without loss of generality we can set $(0) = 0. 
Equation fll5cl) has the following constant solutions 

$ = 0,±1. 

Following the same logic as in Minkowsky case [31,36] since V(<& = +1) = V(<& = — 1) 
we can expect existence of kink-type solutions where scalar field interpolates between 
$ = +1 and $ = —I. 

Iterative procedure (1171) requires initial iteration $ , H Q to be specified explicitly. To 
find such kink solutions it is natural to take 

%(t) = sign(t), H (t) = 0. 

This initial iteration leads to a rapid convergence. It is interesting that our numerical 
experiments show that iterative procedure (fTTj) converges to a kink-type solution if we 
start from any two bounded functions which are positive on positive semi-axis and nega- 
tive on negative one. This observation provides a strong indication for uniqueness of the 
solutions (up to shifts and mirroring). 

The value of £ 2 is determined during level-truncation procedure, its value is 



&y = T~j ~ « 0.96. (24) 



Nevertheless since different values of £ 2 already in Minkowski case result in different 
physical pictures we will systematically consider various values of £ 2 . Note also that the 
case £ 2 = corresponds to p-adic string with p = 3. The value of nonlocal coupling 
constant is fixed by the same procedure it is k = i in our case. 



4 Numerical Results 



4.1 Numerical Solutions for Friedmann Equations 

In what follows for numerical computations we put m p = 1. We used iterative method 
presented in the previous section to solve Friedmann equations for our system. We 
observe the crucial role of the value of £ 2 , as it varies the following physical properties of 
the system are affected (see FigfTJ). 
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For £ 2 = the solution in both $ and H has a monotonic kink shape. This behavior 
is qualitatively similar to the case of Minkowsky metric [31,37,43]. 

There exists a critical value £ 2 SC « 1.18 which determines between types of the $ 
component of the solution - kink or oscillatory with finite period. For £ 2 < £ 2 SC 
the solution has a kink shape with exponentially decreasing oscillations around ±1 
as t — ► ±oo. On the other hand if £ 2 > £ 2 SC the solution converts to oscillations 
with finite periodic] Similar behavior was observed in the case of Minkowsky metric, 
although in that case the critical value was higher « 1.38 [31,37]. 

There is one more critical value £ 2 hapc ~ 0.42 which determines a shape of a Hubble 
function (see FigJ2J). More precisely, as £ 2 grows Hubble function H gets a turning 
point at some positive time t . Moreover, for £ 2 > £ 2 hape H ends up in tending a 
negative (positive) value as t — > oo (t — > — oo). This is a new type of behavior 
which was not present in Minkowsky case. 

Note that the physical value of £ 2 f l24l) is in the region 



f 2 < f 2 < t 2 

S shape >phys S>osC 

Figures [TJ [2] also show how dynamics of the scale factor a(t) 

a(t) = a exp Qf H(r)dr^j (25) 

changes as £ 2 increases. For £ 2 = it has a has a minimum at the perturbative vacuum 
= and increases as <fr tends to nonperturbative vacua ±1. As £ 2 increases the shape 
of a changes and for £ 2 > £ 2 hape the scale factor decreases as (ft — >• ±1 (t — > ±oo). For 
numerical computations we use initial condition a = 1 in 



4.2 Two Profiles of the Hubble Function 

To analyze the late time behavior of Hubble function let us use the so called mechanical 
approximation which was studied for example in [14] and gives realistic qualitative picture 
at least for the late time behavior. Mechanical approximation is constructed by keeping 
in original nonlocal expression only terms with derivatives no higher than second order. 
For (I16bp we get 

H(t)^-— 2 f * dr&dt) 2 - k{d<P) 2 ) = — 2 f dr{k-^){d<j>y. (26) 
m p Jo z m p Jo z 

1 Note that for iterative method which we present in this paper it is essential that solution <I> has 
well defined limits at infinite times. This fact is used in particular for the boundary condition (|20|) . This 
means that this method is not suitable for finding oscillatory solutions with high precision. Oscillatory 
behavior reported in this paper was obtained by vastly enlarging the lattice in t, this did not allow 
finding critical value of £ 2 with more then 3 significant digits. Reliable methods for finding oscillatory 
solutions for the type of equations considered in this paper is an open mathematical problem even for 
Minkowski metric, see [37]. Note though that physically significant solutions in this model turned out 
to be of kink type where iterations rapidly converge. 
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Figure 1: Solutions of the Friedmann equations H and a (left to right) for different 
values of the parameter £ 2 = 0, 0.2, 0.41, 0.43, 0.6, 0.96, 2 (top to bottom). Note that on 
the figures for scale factor a(t) axes cross at the point (t = 0, a = 1). We can see that 
when £ 2 increases the shape of the scale factor changes from parabolic type (a(t) ^ 1 for 
all times) to lump type (a(t) ^ 1 for all times); transition between these shape types is 
illustrated on the Fig. [2j 
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Figure 2: Functions $, H and a (left to right) for £ 2 = 0.41,0.42,0.43 (top to bottom). 
Note how asymptotic behavior changes as we go through £ 2 = 0.42. 

We can see from the expression above that there exists a value of £ 2 which determines the 
change of sign for the Hubble function, more precisely for £ 2 < 2k the Hubble function 
is positive for t —>■ oo (negative as t —>■ — oo), while for £ 2 > 2k the reverse becomes true. 
Thus mechanical approximation gives us the critical value of £ 2 = 2k = 0.25. As already 
described numerical computations give us the value £g hape ~ 0.42 which slightly differs 
from the just found approximate critical £ 2 . This difference manifests the influence of 
higher order nonlocal terms. Nevertheless the analysis above clarifies the phenomenon 
which causes the change of sign for the asymptotic values of the Hubble function. 

4.3 Two Regimes of the Solution 

In this section we will address the mechanism which forms numerically found exponen- 
tially decreasing oscillations and the existence of critical value £ 2 SC . We will also compare 
this behavior with results obtained without gravitational term in action [37]. 
Let us present solution of equation ( 115cl) for the quartic potential V($) = 

(-i 2r D 2 n + l)e 2fcI) H$ = $3 (27) 

as a sum 

$(t) = $o(t)+x(t), (28) 

where $o denotes solution of equation fll5cl) for £ 2 = 0. Substituting (|28|) to (I15cl) and 
leaving only linear terms in x we get 

(-£ 2 D H + l)e 2fcD H($ + x ) = $3 + 3 $2 X . (29 ) 
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Using the fact that $o is the solution of equation (I15cp for £ 2 = we can write equation 
for x(t) which in linear approximation has the form 

{-e^l + l)e 2k ^X = 3$ lx + e%e 2kv2 ^ . (30) 

Using the fact that operator e 2fcI) H acts as an identity on constants we can write large t 
approximation 

(-e% + l)e 2k ^X = 3X, (31) 

were we used asymptotical properties of solution $ - 

To carry out harmonic analysis for equation (13T1) we will consider eigenfunctions of the 
D'Alembertian operator D^vx = —\v\. Using that e 2kT ^v\ = e~ 2kX v\ and expanding x 
in v\ equation fl3Tj) leads us to 

(£ 2 A + l)e- 2fe \; A = 3t; A . (32) 

Our goal is to analyze which values of £ 2 allow real-valued A, i.e. non-damping oscilla- 
tions take place. Considering equation (1321 as an equation for complex variable A with 
parameter £ 2 we obtain that there is a minimum value of £q ~ 1-77 for which A is real. 
So for £ 2 < £ 2 we have solutions with nonzero imaginary parts which result in vanishing 
solutions while £ 2 > £q leads to oscillatory regime. Equation ( 132]) has exactly the same 
form in case of Minkowsky space which was considered earlier [37] (functions v\ are dif- 
ferent as they depend on H) . As we can see method discussed in this section provides us 
only qualitative explanation of the changing of the regimes, the value of critical £ 2 found 
here only approximately reproduces numerically obtained £ 2 SC ~ 1.18. It is interesting 
to note that Minkowski case leads to higher value of critical £ 2 [37], i.e. it appears that 
"friction" H(t) does not damp oscillations but au contraire increases them. 



5 Cosmological properties 

It is interesting to note, that initially we considered a model with potential V((p) = 
— \(f) 2 + + see Fig. [3^). Nonlocal term contributes to kinetic term, most importantly 




Figure 3: W and —W shaped potentials. 



it changes its sign [11], or equivalently changes the sign of potential (kinetic term is 
positive in this case), see Figj3l In this context it is interesting to study already mentioned 
mechanical approximation which can illustrate such an effect. 
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Solutions of (fl5j) lead to interesting cosmological properties. Let us consider them 
for physically interesting values of £ 2 and see what kind of physical behavior they result 
in. We would like to study dynamics of the following physical quantities. The state 
parameter which is defined as usual 

n 

(33) 



w 



P 



P 



or in the terms of Hubble parameter 

» = -i-!£. <*) 

will be presented on the figures below along with deceleration parameter defined as 

(35) 



. . aa 

q(t) = — • 

a 



5.1 Effective mechanical potential 

In order to better understand which behavior of the system we might expect let us first 
study equation (jUJ) in local approximations. 

If we ignore the nonlocal operator e 2fc ^H in © we get 



(-£ 2 Dj?j + 1)$ = $ 3 . 
Let us reproduce action which leads to the equation above 



S 



m £ 2 1 1 

2 2 9 2 4 



This action describes a particle moving in the system with friction H(t) 

£ 2 <9 2 $ = $ - $ 3 - 3#£ 2 d$, 

H = -il(a$) 2 . 



(36) 

(37) 

(38a) 
(38b) 




a . 



b). 



Figure 4: Anharmonic oscillator with W-shape potential; a) periodic scalar field trajec- 
tory located around $ = 1; b) the solution for the Hubble function H(t), £ 2 = 0.4. 
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Figure 5: Scalar field $(£) (upper row) and Hubble function H(i) (lower row) for the 
system (1411) for different values for £ 2 parameter £ 2 = 0, 0.2, 0.4, 1 from left to right 
correspondingly. 



As we can see here Hubble function is explicitly negative and thus we obtained a system 
which is unusual from mechanical point of view - it is an anharmonic oscillator with 
negative friction (see Fig. H]). Such a behavior looks similar with the behavior for nonlocal 
systems studied in [14, 15] in the Minkowski background, when the trajectories of the 
scalar field exceed the extremum of the potential, meanwhile the energy is conserved. 

Now let us consider mechanical approximation by accounting in the equation (Q only 
derivatives no higher that second order, we get 



((2k - f)d 2 + 1)$ = $ 3 - 3(2A; - f)H(t)d<5>. 
Equation (1391) could be obtain from the following action 

which leads to the following local Friedmann equations 

$ _ $ 3 



(39) 



(40) 



<9 2 $ 



H 



(2k - a 
(2k - a 

ml 



- 3Hd$, 



(d$f 



(41a) 
(41b) 



This system of equations describes particle moving in the following potential 



2$ 2 - $ 4 
A(2k-i 2 )' 



As we can see for £ 2 < 2k flip of the potential takes place. In particular, for £ 2 = we 
get an equation describing particle moving in the — VT-shape potential (see FigfSb). For 
£ 2 > 2k the potential has VT-shape (see Figj3^i). This simple remark leads to important 
physical consequences, we might expect two different regimes of the solution. Solutions 
for scalar field and Hubble function for different values of £ 2 are presented on Figj5j We 
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see two regimes - decreasing oscillations with positive Hubble function and increasing 
oscillations with positive Hubble function. 

It is important to note that while mechanical approximation helps us to illustrate or 
explain qualitatively some physical consequences, it is a rather rough approximation and 
does not necessarily capture some nonlocal effects. This situation is already depicted 
by the fact that in full nonlocal model and its mechanical approximation the change of 
regime of the solution corresponding to W and — W shapes of the potential happens for 
different values of £ 2 , they are of the same order though. 

5.2 Rolling tachyon in the FRW background and the effective 
cosmological constant A' 

Let us study the value of effective cosmological constant A'. Formally A' enters action 
as a correction of D3-brane tension T. For rolling solution in Minkowski background 
it is known [11,31,44] that the value of T is essential for calculating stress tensor, but 
as it does not enter equations of motion it is not essential for the existence of solution. 
Moreover according to Sen's conjecture [11,31,44] 

^Minkowski = — V(<& = ±1) = -. 

Such value of T in Minkowski background corresponds to zero cosmological constant in 
the non-perturbative vacuum $ = ±1. 

In case of FRW metric we represent D3-brane tension as T = T m i n kowski + A' where A' 
is determined uniquely for each field configuration (situation is the same in local theories). 
Note that A' does not enter equations f|T6|) . in fact we determine its value from f)15al) . 
Different values of A' for different £ are presented on table [TJ Note that we worked with 
dimensionless variables. If we return to initial notations and take into account physical 
constants in action ([6]) we can get realistic value of A because generally speaking the 
string scale does not coincide with the Plank mass [11, 17, 19]. 



e 


A' 





0.036938 


0.2 


0.010379 


0.41 


0.000022 


0.42 


0. 


0.43 


0.000023 


0.6 


0.007766 


0.96 


0.077848 



Table 1: Values of A' for different £ 2 . 
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Figure 6: The state and deceleration parameters for the case of p-adic cosmological model 

(e = o). 



5.3 P-adic cosmological model, £ 2 = 

Recently p-adic string models started to attract a lot of attention being used as toy 
cosmological models. In [19] p-adic string model [31] was considered in the context of 
inflation and approximate solutions of the fully nonlocal p-adic string theory coupled to 
gravity were constructed. In section 13.11 we presented a new algorithm for construction 
of precise solutions for these models. While iterative procedure of section 13.11 converges 
for any odd p (for rigorous mathematical proof for the simpler Minkowski case see [43]) 
here we will analyze cosmological properties of the solution for the case p = 3. This 
case is specifically interested since as mentioned earlier it corresponds to zero mass of the 
tachyon field, £ 2 = 0, in level truncated fermionic string field theory (SFT). 

We can see on Fig{T]and Figj6]that scalar field and background solutions are monotonic 
functions which tend to constants as t — > ±oo. Since these solutions do not have turning 
points (oscillations) there is no crossing of the cosmological state parameter barrier (w = 
— 1), it is negative and approaches to —1 from below (Figj6]). For these solutions we have 
nonsingular accelerating Universe with a bounce. As an illustration we can note that 
behavior of the scale factor obtained could be well approximated by a a pprox(£) = e Hot for 
large t > and by a approx (t) = e~ Hat for large t < 0, where H = lim^oo H(t). Thus we 
get de Sitter type solution for flat FRW background for t — ► +oo and anti de Sitter as 
t — > — oo. It might be useful to remark that bouncing cosmology is a subject of of recent 
investigation which include the Ekpyrotic [34], Pre-Big-Bang [45] and higher derivative 
modification of Einstein gravity scenarios [46] . 

Deceleration parameter q is negative and approaches —1 from below as t — > ±oo. 



5.4 SFT cosmological model, £ 2 = £p hys 

Today SFT is the strongest candidate for real physical theory of nature. As we mentioned 
earlier level truncation in non-sypersymmetric fermionic SFT leads to action fl7|) with po- 
tential (1211) with £ 2 = £p hys (1241) [11,31,33]. The algorithm presented in section [37TI allowed 
us to study this model coupled to FRW backgroud without explicit approximations. 

As shown on Figfj] and FigJT] in this case scalar field solution has turning points 
(decreasing oscillations) and tends to stationary solutions ±1 as t — » ±oo. The shape of 
Hubble function is essentially different as compared to the case £ 2 = 0, indeed it has a 
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Figure 7: The state parameter, its fine structure and deceleration parameter for the case 
of SFT cosmological model (£ 2 = £p hys ). 

clear maximum and tends to negative constant for t > 0, lim^-too H(t) = =f0.15. Such 
behavior is in agreement with analysis presented in the previous section since £ph y s > Cose- 
The state parameter w(t) tends to — 1 as t — > ±oo and has exponentially decreasing 
oscillations around w — — 1 barrier. Note that models which allow for crossing of w — — 1 
barrier are a subject of many recent investigations [17,47-53]. Here we observe this 
behavior in nonlocal model with only one scalar field as in [17]. We also see that we 
have a decelerating phase at late times, this property was not predicted by approximate 
methods [11]. 

6 Conclusions 

In this paper we have studied the dynamics of nonlocal cosmological models driven by 
String Field Theory. These models have an infinite number of derivatives and are charac- 
terized by positive constants k, £ and A' which determine the shape of the solution. We 
have developed a new method for solving nonlocal Friedmann equations. Such equations 
contain infinite number of derivatives and form a new class of equations in mathematical 
physics which recently started to be discussed in literature [15,18,19,37,41-43,54]. 

To study cosmological properties we have investigated the behavior of nonlocal models 
in the FRW background. We are especially interested in two physical cases: p-adic 
cosmological models and level-truncated SFT model. In p-adic cosmological model we 
obtained nonsingular bouncing solution with w parameter approaching the —1 barrier 
from below. SFT model has nonsingular solution for which w parameter crosses — 1 
barrier. We also discussed the possibility of obtaining realistic cosmological constant 
from considered above nonlocal cosmological models. 
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